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c-Axis Penetration Depth in the Cuprates: 
Additional Evidence for Incoherent Hopping 

R. J. Radtke and V. N. KosturO 

Center for Superconductivity Research, Department of Physics, 
University of Maryland, College Park, Maryland, 20742-4111 
(E-print cond-mat/9509nnn, received February 1, 2008) 

Measurements of the c-axis penetration depth A c in the cuprates reveal a low-temperature T 
dependence which is inconsistent with simple models of coupling between the Cu02 layers. In this 
paper, we examine whether a model based on incoherent hopping between the layers can account for 
this low-T behavior. We compute A c directly from linear response theory and compare our results 
with recent experimental measurements on YBa2U307_<s as a function of temperature and doping. 
We find that the data can be reproduced within this model providing the inter-layer scattering 
is anisotropic and the pairing is d- wave. In addition, our calculations demonstrate that is 
proportional to the c-axis critical current, which seems to be a generic feature in weakly coupled 
layered superconductors. 



Introduction. Studies of the c-axis properties of the 
cuprate superconductors have revealed many unusual fea- 
tures in both the normal and superconducting states. ||] 
In the normal state, the c-axis resistivity p c (T) can show 
either a metallic or a semiconducting temperature de- 
pendence as a function of the doping, while the in-plane 
resistivity is always metallic, ffl In the superconducting 
state, strong evidence exists in a variety of cuprates that 
the CuC>2 planes are coupled by the Josephson effect, 
in contrast to the behavior of conventional superconduc- 
tors. Additionally, c-axis penetration depth mea- 
surements on YBCI2CU3O7-5 (YBCO) crystals have be- 
come available |5|,^j which agree with neither the conven- 
tional theory nor with each other. 

Recent work has been directed toward accounting for 
these anomalous features within a theoretical framework 
based on incoherent quasiparticle hopping between the 
CuC>2 planes. Q In this picture, the coupling between 
nearest-neighbor CuC>2 planes is so weak that quasipar- 
ticle transport occurs through uncorrclatcd tunneling (or 
hopping) of the quasiparticles. The origin of this incoher- 
ence is unclear, but it may be related to strong intra-layer 
scattering, a non- Fermi-liquid ground state within the 
layers, H dynamic inter-layer scattering, [Bj or strong 
electromagnetic fluctuations. Jl0| Based on these ideas, 
several phenomenological models have been developed to 
study the implications of this incoherence 



I! 



II and 



of 



have had some success in describing both p c {T) 
and the c-axis critical current j c . [ fTsjfhj 

In this paper, we extend the calculations 
Rcfs. 12 14 by deriving expressions for the electromag- 
netic penetration depth along the c-axis directly from 
linear response theory and compare our results with re- 
cent experiments. || Following Ref. |Q, we postulate the 
existence of three mechanisms for quasiparticle transport 
along the c-direction: a direct hopping induced by over- 
lap of the quasiparticle wave functions, impurity-assisted 



hopping due to disorder in the intercalating layers sepa- 
rating the CuC>2 planes, and boson-assisted hopping due 
to, for example, certain optical phonons in the cuprates. 
We find that the inverse c-axis penetration depth 1/A C 
is the sum in quadrature of 1/A C due to each mecha- 
nism in isolation, and we derive expressions for these 
terms. Our calculations demonstrate that 1/A;: oc j c as 
hypothesized earlier |l3],[l4| based on calculations within 
the Lawrence-Doniach model, |1| indicating that this 
Josephson-like relation between j c and A c is a general 
property of weakly coupled layered superconductors. Fi- 
nally, we examine the temperature dependence of the 
contributions of each interlayer hopping process to the 
penetration depth for both s- and d-wave superconduc- 
tors. Comparing with the data of Bonn et al, || we find 
that the doping and temperature dependence of the ex- 
perimental A^(0)/A^(T) is qualitatively consistent with 
d-wave pairing and anisotropic interlayer scattering. 

Penetration depth in the incoherent hopping model. To 
begin, consider electrons on a three-dimensional tetrago- 
nal lattice governed by the Hamiltonian H c \ — J2 m + 
H±. The lattice is defined by the primitive vectors ax, 
ay, and di with d > a, so that it resembles a stack of 
layers, which we index by m. H m is the Hamiltonian for 
the motion in the mth layer and H± describes the inter- 
layer coupling. For the purposes of this paper, it is not 
necessary to specify H m completely; we need only assume 
that it becomes a BCS-type superconductor below a criti- 
cal temperature T c . The interlayer coupling Hamiltonian 
may then be written in terms of the annihilation and cre- 
ation operators for the quasiparticles obtained from H m 
at site i = (i x , i y ) in layer m with spin projection a, Ci ma 
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The terms in this Hamiltonian represent interlayer hop- 
ping due to overlap of the quasiparticle wave functions 
(parameterized by t±), impurity scattering (modeled by 
the random variable Vi m ), and bosonic scattering (writ- 
ten in terms of the bosonic field operator <p im and its 
coupling strength to the quasiparticles <7j TO ). Jl6| 

As discussed in the Introduction, a significant number 
of experiments suggest that the interlayer coupling in the 
cuprates is in general weak and in particular incoherent. 
A fully microscopic theory of this incoherence is beyond 
the scope of this article, but its effects can be simulated 
by performing calculations to second order in H±. This 
approach is similar in spirit to that of Ioffe et al. |l7| ] 
and has been used successfully to study the normal-state 
resistivity Jl^,[l3| and the c-axis critical current. [ p~3|Jl4| 
We take this approach in what follows and refer the inter- 
ested reader to Ref. [0 for a more detailed description of 
the perturbation theory and the assumptions underlying 
it. 

To compute A c within this incoherent hopping model, 
we examine the current response to a weak vector po- 
tential Ai m (t)i in linear respone theory. As described by 
Peierls, |l8| the presence of the vector potential induces a 
phase e l ( ed / hc ) A im j n to Eq. (|l|), where d is the inter-layer 
spacing, e is the electronic charge, and c is the speed of 
light. For weak fields, we may expand the exponential to 
second order in Ai m and thereby obtain an expression for 
the total charge current density: jl9) j 
where the paramagnetic current 
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proportional to the Fourier transform of a retarded 
current-current correlation function and 

e 2 d 1 



D 



Q,Q' 



% 2 a 2 N 



-i(Q-Q')-Ri 



■%m+l,cr L ''imtr 



(7) 



the diamagnetic term. In these expressions, N is the total 
number of sites in the lattice, j'q P = J2im e ~ lC> Rim Jim j 
Q = (l, Qz), and R^ m = (r^, md) is the position vector of 
the lattice site indexed by im. (Throughout this paper, 
we use capital letters to denote 3D vectors and small 
letters to denote 2D vectors) . Note that the conductivity 
is not diagonal in wave vector due to the impurity term; 
after impurity averaging, translational invariance will be 
restored and we will find <7q q, (lu) oc <5q,q'- 

To proceed further we must account for the incoherent 
nature of charge transport along the c-axis by expand- 
ing both IIqq, and Dqq, to second order in the inter- 
layer hopping amplitude ti m . Thus, the paramagnetic 
current-current correlation function, evaluated using the 
standard diagrammatic techniques in Matsubara space, is 
represented by a bare particle-hole bubble (vertex correc- 
tions being of higher order in ti m ) with purely infra-layer 
Green's functions renormalized only by the infra-layer 
self energy. Although the diamagnetic term does not 
play any substantial role in the optical conductivity be- 
cause it is independent of frequency, it is important for 
the penetration depth calculation and must be treated 
carefully. In fact, evaluating this term perturbatively 
yields a form nearly identical to Uq q, except for fac- 
tors at the vertices, as we shall see below. Working in 
Matsubara space with Nambu Green's functions for the 
intra-layer propagators ]2l[ | and assuming that the layers 
to are identical, we obtain 



and the diamagnetic current 
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In the usual way, [£0| we expand the expectation value 
of the total current to linear order in Aj m , Fourier trans- 
form the result in space and time, and use Maxwell's 
equations to relate the electric field Ei m (t)z to the 
vector potential. This procedure gives jQ Tot (^) — 
Eq' ^q.q'H^Q'H* where 
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where Nu is the number of lattice site in a single layer, 
f 3 is the third Pauli matrix, and i q = ^2 i e~ l<i ' ri ti m is 
assumed to be independent of to. (kg = 1 throughout 
this paper, and the rest of the notation is standard.) 

The correlation function (r r f q (r)t q /) represents a com- 
bination of an impurity average over Vi m and a thermo- 
dynamic average over the bosonic field <f>i m . Since we go 
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to second order in the hopping amplitudes and take the 
mean of Vi m to be zero, this correlation function decom- 
poses into a sum of three components: a direct hopping 
term, an impurity-assisted hopping term, and a boson- 
assisted hopping term. The conductivity is therefore a 
sum of three terms, as obtained previously. Since the 
penetration depth is given by 
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after impurity averaging, this decomposition allows us to 
write 



Total 
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and consider each term separately. 

Direct contribution. For the direct term, ti m 
Eq. (||)] , from which 
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The resulting conductivity is, as expected, diagonal in 
wave vector. Inserting Eqs. (jl^)-( 14) into Eq. (|l0|), we 
have 



direct 
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where Fk(iwj) is the Gor'kov propagator. pl| 

Several features of this result are noteworthy. First and 
foremost, Eq. ( |l5| ) has the same form as the expression 
for the c-axis critical current; specifically, (1/A;:)direct 
jdirect^ |TJjj ^ g we ghoji see ^ ^his conclusion holds for the 
other two processes as well. This behavior arises from 
the cancellation of the normal (fo and T3) components of 
the quasiparticle propagators between the paramagnetic 
and diamagnetic terms. This cancellation is a special 
property of the c-axis charge transport in weakly cou- 
pled layered superconductors and reproduces the results 



obtained by Bulaevskii and Clem based on the Lawrence- 
Doniach model. |l5| 

Second, if we take the BCS form for the Gor'kov prop- 
agators -Fk(iwz) = Ak/((iw;) 2 — -^k)> m l nc limit where 
the superconducting energy scales are smaller than the 
intra-layer electronic energy scales, we obtain 
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where iV(0) is two-dimensional density of states at the 
Fermi surface (for a parabolic spectrum ek = Ti 2 k 2 /2m — 
Ep, N(0) — m/2nh 2 ), Ak is the superconducting gap 



function, E* = \fe^ -r ^ k 



A^, and the angle brackets de- 



note a normalized Fermi surface average over the indi- 
cated variable. Comparing with the standard BCS re- 
sult, [§2j we see that this term has the same temperature 
dependence but a different magnitude, (lj] This behav- 
ior arises from a combination of the dimensionality of 
the layers and the conservation of the transmitted wave 
vector. (H| 

Impurity-assisted contribution. The impurity-assisted 
term is computed in an analogous way. In this case, 
Urn — * Vi m [cf- Eq. (pj)] and the impurity average yields 
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The penetration depth follows from these relations as 
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Again, the penetration depth has the same form as the 
corresponding critical current calculation: (1/A 2 ) im p oc 
j™ p . @ In the BCS limit taken above, Eq. @ may be 
written 

32n 3 e 2 v 
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where v = a 2 d is the unit cell volume. 

The contribution of the impurity-assisted component 
to the total penetration depth depends on both the pair- 
ing symmetry and the anisotropy of the matrix element. 
If the pairing is isotropic (Ak = A), this expression re- 
duces to the Ambegaokar-Baratoff result for SIS tunnel 
junctions. p3| However, if the pairing is d-wave (or any 
pairing state with (Ak) k = 0) and the matrix element 
is isotropic (V 2 k-k' = V 2 ), (V^c)imp vanishes by sym- 
metry. 0,H This result is a special consequence of an 
isotropic scattering matrix element and does not hold 
when anisotropy, which should be present in any real 
material, is included. 

To see the effects of anisotropy on the impurity-assisted 
contribution to the c-axis penetration depth, consider a 
scattering matrix element of the form 
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When 5k — » oo, V r2 k-k / becomes isotropic, and we re- 
cover the results discussed above. Physically, this limit 
corresponds to a total randomization of the wave vector 
of the scattered quasiparticle and is analogous to diffuse 
transmission through a tunnel junction. In the opposite 
limit where 5k/ kp — > 0, the scattering matrix element 
preserves the direction of the scattered quasiparticle and 
is analogous to specular transmission through a tunnel 
junction. In this limit, Eq. ( ^l|) reduces to a generaliza- 
tion of the Ambegaokar-Baratoff |E3] result: 



Ag(0) 
A 2 (T) 



<5fc^0 
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where Ao = (|At(T = 0)|) k is the average of the absolute 
value of the gap function over the Fermi surface at zero 
temperature. This form yields a finite contribution to 
the penetration depth for both s- and <i-wave pairing due 
to the wave- vector-conserving nature of the scattering. 

For intermediate 5k/kp, we can evaluate (l/A 2 ) imp 
from Eqs. (^i"|)-(p2|) for a d-wave superconductor with 
the results shown in Fig. [j]. At zero temperature, we see 
from the inset to this figure that the magnitude of the 
penetration depth depends strongly on 5k/kp, in keeping 
with the fact that (Aj 2 (0))i mp vanishes in the 5k — > oo 
limit. As seen from the main figure, however, the tem- 
perature dependence of the penetration depth ratio is not 
strongly affected by the amount of scattering anisotropy: 
the penetration depth ratios for 5k/kp = 0.01 and 2.0 
lie very nearly on top of each other. Moreover, we 
see that the temperature dependence of the anisotropic 
impurity-assisted scattering is distinguishable from both 
the conventional BCS result and the s-wave Ambegaokar- 
Baratoff result. 

Our calculations therefore suggest two experimen- 
tal signatures of anisotropic impurity-assisted interlayer 
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FIG. 1. Normalized contribution to the c-axis pene- 
tration depth from impurity-assisted interlayer hopping 
(Ac(0) / \c(T)j . as a function of the reduced tempera- 
ture T/T c in a d-wave superconductor with slight (Sk/kF = 
0.01, solid line) and significant (Sk/kF = 2.0, dotted line) 
non-conservation of wave vector [cf. Eq. (0)]. For compar- 
ison, the d-wave BCS (dashed line) and the s-wave Ambe- 
gaokar-Baratoff results (dot-dashed line) are also shown. 
Inset: Square of the c-axis penetration depth at zero temper- 
ature Xc(0) as a function of the spread in transmitted wave 
vector 5k jkp [cf. Eq. (0)] normalized to its value at Sk/kF 
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scattering: (1) the magnitude of the c-axis penetration 
depth should depend strongly on the type and quantity 
of disorder, but the temperature dependence should be 
largely insensitive to these factors, and (2) the tempera- 
ture dependence of the penetration depth ratio should be 
between the BCS and Ambegaokar-Baratoff predictions. 

Boson- assisted contribution. The boson-assisted hop- 
ping term is more involved due to the inelastic scat- 
tering (i.e., the retardation effects) induced by the bo- 
son. The hopping amplitude t irn — > £\ g l -^ m 4>j m 
[cf. Eq. (0)], and so the thermodynamic average re- 
turns a factor of the boson propagator Dq m (r) — 



(T T [4>v 



relation function in Eqs. (pf) 
(T T £ q (r)t q . 



I incl 



NnS. 



1,-1 



i(0)]), allowing the cor- 
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|g qm | 2 ^ qm (r). (24) 



If the layers are identical, we may drop the index m and 
compute the conductivity and the penetration depth as 
in the preceding two cases, giving 
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As before, this expression is proportional to the corre- 
sponding one for the critical current. Jl4|] In the BCS 
limit, we can introduce the spectral representation of the 
boson propagator familiar from the theory of supercon- 
ductivity, Hi 

B k , k -(fi) =-- | 5k -k<| 2 ImD k _ k ,(ft), (28) 
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and write the result as 
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For a simple Einstein phonon with a structureless cou- 
pling constant (-B k ,k'(^) & ^(^ — ^o))j we nn d that 



resembles the Ambegaokar-Baratoff form for 



s-wave pairing and Slo > Ao , but vanishes when the pair- 
ing is ci- wave. 14 As with impurity-assisted scattering, 
however, anisotropy in -B k ,k' (f2) will in general be present 
and will have a profound effect on (l/X 2 ). . in d-wave 
superconductors. Thus, we need to quantify the effects 
of anisotropy in this case. 

The anisotropy in the boson-assisted scattering may 
occur due to wave-vector dependence in either the 
electron-boson coupling constant g q or the dispersion of 
the boson itself, so it is difficult to obtain a simple form 
for -B kjk '(f2) similar to Eq. (p2|). We therefore take an 
alternate approach and expand the wave vector depen- 
dence of all quantities in Eq. ( |29| ) in terms of Fermi sur- 
face harmonics .Fz,(k). pj| We further simplify our re- 
sults by taking an Einstein spectrum for each component 
of B, B L>L ,(Q) = (Cl L , L >/2N(0))6(Q-n LtL ,). This 
procedure gives (1/A^) ine i = J2l //(V^cKi/ with 
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and At{iuJi) = ^F L (k)A k / v / wf+Ajy . For concrete- 

ness, we take a cylindrical Fermi surface within each 
layer, which gives ■Fi(k) = \[2 cos(L<fi) , where k = 
hp (cos <p, sin <j>). 

We can make several general statements by looking at 
the form of Eq. (^0|). First, although the total (l/A^inei 
must be positive, the individual components (1/\ 2 )l,l' 
may have either sign because of the difference in the 
signs of Al and A^i . Second, for isotropic s-wave pairing 
(A k = A), only the L = L' = component is non-zero, 
and we recover the results of Ref. iQ after analytic con- 
tinuation. For d-wave pairing, on the other hand, only 
components with (L,L') = (An + 2, An' + 2) (n and n' 
integers) are non-zero. 

These features are exhibited in Fig. ||, which shows 
(1/\ 2 )l,l> for both s- and d-wave pairing computed nu- 
merically from Eq. (^0|). In these calculations, we take 
Hl.L' = Qph — 40 meV to reflect the expected impor- 
tance of this phonon mode in c-axis transport. p6| For 
isotropic s-wave pairing, only the (L, L') = (0, 0) com- 
ponent contributes, and we obtain a result which is sim- 
ilar to, though slightly different from, the Ambegaokar- 
Baratoff result as long as the boson involved has an en- 
ergy scale larger than the maximum of the gap function. 
fll4f For d-wave pairing, many (L, L') components are 
non-zero and some are negative, as suggested above. 

Examining the d-wave curves more closely reveals sev- 
eral interesting features. At low temperatures, the non- 
zero components do not have the linear-in-T behavior 
typical of d-wave superconductors, but seem to follow 
a higher power law. This result is consistent with the 
fact that the nodes in a d-wave gap function must give 
rise to a power law in the penetration depth, but we 
see that the exact exponent for this power law is deter- 
mined by the mechanism responsible for the hopping and 
the symmetry of its matrix elements. At higher tem- 
peratures, the contribution from higher harmonics (L, 
V > 2) vanishes. This feature results from the orthog- 
onality of the Fermi surface harmonics in the following 
way. As T — > T c , the magnitude of Ak decreases and 
AL(ioj n ) ~ (i*z,(k)Ak) k /\u n \- Because the gap function 
is pure <i-wave (L = 2), A k cx ^(k) and so AL(ito n ), 
and hence (l/X 2 ) L L ,, goes to zero when L, L' > 2 much 
faster than for the L = V = 2 component. We have also 
computed the magnitudes of the penetration depth com- 
ponents (l/X 2 ) LL , and find that they depend strongly 
on (L, L'), becoming considerably smaller with larger L 
and L' . Thus, one expects that the (L, L') — (2, 2) com- 
ponent will dominate the boson-assisted contribution to 
the c-axis penetration depth in a d-wave superconductor. 

The conclusion one can draw from this analysis is that 
boson-assisted hopping in a d-wave superconductor can 
contribute to the c-axis penetration depth if the coupling 
is anisotropic or if the boson has some dispersion. More- 
over, the resulting contribution to the penetration depth 
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FIG. 2. Normalized matrix elements of the c-axis pen- 
etration depth in the basis of Fermi surface harmonics 

(A^ 2 (T)) / (A ( T 2 (0)) i as a function of the reduced 

temperature T/T c for a d-wave superconductor and (L, I/) 
= (2,2) (solid line); (6,6) (dashed line), (2,10) (dot-dashed 
line); (10,10) (long dashed line); (6,10) (dot-dot-dot-dashed 
line); and (2,6) (dotted line). Due to the symmetry 
of the the order parameter, only matrix elements with 
(L,L') — (4n + 2,4n' + 2), n and n' positive integers, are 
non-zero. For comparison, the s-wave, (L,L') = (0,0) result 
is plotted as the thick solid line. See text for details. 



is not necessarily linear in T at low T, contrary to what 
one expects and similar to what is observed experimen- 
tally. 

Relation to experiment. In order to connect our re- 
sults more firmly to experiment, we compare some rep- 
resentative curves in our model to the experimental A c 
data of Bonn et al. @ on YBCO in Fig. |. The dashed 
line shows the BCS prediction for a d-wave superconduc- 
tor: it is clearly inconsistent with these data. The solid 
and dotted curves show the penetration depth ratio for 
purely disorder-mediated hopping and a linear combina- 
tion of disorder- and boson-assisted hopping. We see that 
the penetration depth ratio computed within the incoher- 
ent hopping model is qualitatively consistent with these 
data: despite having a c£-wave order parameter, the low- 
temperature behavior shows a much smaller slope, and 
the penetration depth ratio falls much more slowly than 
BCS theory predicts. 

We are also able to understand the systematic change 
in the temperature dependence of the penetration depth 
ratio with de-oxygenation. Several experimental and 
structural features of YBCO, along with an analy- 
sis of its resistivity within our incoherent hopping model, 
| p2| pT| suggest that the fully oxygenated compound is 
both "cleaner" and more three-dimensional than the de- 
oxygenated samples. In particular, we find that the 
boson-assisted hopping is more pronounced in the de- 
oxygenated compound and is, in fact, responsible for 
the upturn in the c-axis resistivity with decreasing tern- 
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FIG. 3. Normalized c-axis penetration depth ratio 
A 2 (0)/A 2 (T) as a function of the reduced temperature T/T c 
for a d-wave superconductor with coherent (BCS-like) in- 
ter-layer coupling (dashed line), incoherent coupling due to 
disorder-assisted hopping with Sk/kF = 0.01 (solid line; cf. 
Eq. ([22|)), and incoherent hopping from a linear combination 
of 80 % disorder- (Sk/k F = 0.01) with 20 % boson-assisted 
hopping (dotted line; cf. Eq. (||) with (L, L') = (2,2)). Also 
shown are the data of Bonn et al. || for T c = 93 K (open 
squares) and T c = 60 K (open triangles) YBa2Cu307_<s . Ob- 
serve that effect of adding boson-assisted hopping is qualita- 
tively the same as de-oxygenation in the experimental data: 
the penetration depth ratio is larger at larger temperatures. 



perature. |12|-|14| We may therefore expect that the o 
axis penetration depth ratio will have a larger contri- 
bution due to the boson-assisted component in the de- 
oxygenated samples while the disorder-assisted compo- 
nent remains relatively constant. The net result is that 
the penetration depth ratio will fall more slowly with 
temperature in the de-oxygenated samples, and this is 
what is observed experimentally. 

Summary. Viewing the Cu02 layers in the cuprate su- 
perconductors as incoherently coupled explains the fea- 
tures of the anisotropic resistivity in these materials 
|]l2| , ^3| and makes definite predictions for the c-axis criti- 
cal current j c . Q] We have extended this theory by pro- 
viding a microscopic derivation of the c-axis penetration 
depth A c and find that j c oc 1/A^, which seems to be a 
generic property of superconductors consisting of weakly 
coupled layers. We have also examined the contri- 
butions to A c arising from direct, impurity-assisted, and 
boson-assisted hopping. We find that anisotropy in the 
scattering matrix elements or boson dispersion are re- 
quired in order to obtain a finite penetration depth in d- 
wave superconductors in the absence of direct hopping, 
although no such restriction arises for s-wave pairing. 
By computing the temperature dependence of the pene- 
tration depth ratio X^(0) / X^.(T) for d-wave pairing and 
anisotropic scattering, we find that the experimental c- 
axis penetration depth ratio as a function of both tem- 
perature and doping is qualitatively consistent with this 



G 



model for at least some measurements of A c in YBCO. 
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